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1.0 INTRODUCTION 


The following is the final report on the comparison of 
analytical and experimental pressure measurements on four oscil- 
lating helicopter rotor blade tips. Each blade tip is considered 
as a semispan wing oscillating in pitch about the quarter-chord 
axis. The four blade-tip planforms investigated include a swept, 
a rectangular, a tapered and an ogee shape. Figure 1. 

Two computer programs, VSAERO-TS and VSAERO-H, were de- 
veloped for computing the unsteady subsonic aerodynamic charac- 
teristics of arbitrarily shaped wings oscillating in pitch, bgth 
programs have a common basis (Figure 2) in program VSAERO, which 
is under continued development (Refs. 1 through 4) for the analy- 
sis of steady, non-linear aerodynamic characteristics of arbi- 
trary configurations. Program VSAERO-TS is a time-stepping 
analysis capable of treating large amplitude motions while pro- 
gram VSAERO-H uses harmonic wake and small amplitude assumptions. 

The basis of the computer program is a surface singularity 
panel method using quadrilateral panels on which doublet and 
source singularities are distributed in a piecewise constant 
form. The panel source values are directly determined by the 
external Neumann boundary condition controlling the normal com- 
ponent of the local resultant flow: the doublet values are 
solved after imposing the internal Dirichlet boundary condition 
of zero perturbation potential at the centers (underside) of all 
the panels simultaneously. Surface perturbation velocities are 
obtained from the gradient of the doublet solution, while field 
velocities are obtained by direct summation of all singularity 
panel contributions. The details of the mathematical formulation 
are presented in a separate theory document (5). 

The program is written in standard FORTRAN IV and has been 
developed on the CDC Cyber 175 computer. Minor changes to the 
code allow it to run on an IBM 3033 computer. The time-stepping 
version of the unsteady code allows up to 1,000 panels (i.e., 
unknowns), while VSAERO-H, which uses complex variables, allows 
up to 240 panels. The core requirement (octal) for either pro- 
gram is 234 K. 

The experimental results were obtained by the Institut fur 
Aeroelastik in the DFVLR 3 x 3 m tunnel in Gottingen, West 
Germany. The work was performed under an international coopera- 
tive agreement between the DFVLR and NASA. 

In the following sections, the correlation between the com- 
puted and the experimental chordwise pressure distribution for 
the four tip shapes is presented. For the swept tip case, a 
complete set of cases covering the range of 0.1, 0.2 and 0.3 
reduced frequencies and angles of attack of 0, 4, 8 and 12 de- 
grees are included. Only a limited number of cases are included 
for the other three tip cases. The scope of the present correla- 
tion excludes the presence of extensive dynamic separation 
effects. 
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Figure 2. 


VSAERO Unsteady Program Development. 






2.0 SWEPT TIP 


VSAERO-TS is used to compute several cases for swept tip 
wingsi The following cases of comparison between the computed 
and the DFVLR experimental pressure distributions are presented. 


Figure No. 

a o 

(degrees) 

a i 

(degrees) 

a) 

3 

4 

0.710 

0.1 

4 

12 

0.704 

0.1 

5 

4 

0.707 

0.2 

6 

8 

0.707 

0.2 

7 

0 

0.714 

0.3 

8 

4 

0.714 

0.3 

9 

8 

0.717 

0.3 

10 

12 

0.704 

0.3 


Table 1. Comparison between Computed and DFVLR Experimental 
Pressure Distributions for the Swept Tip Wing. 


The cases included in this section cover the angle of attack 
(mean) range of 0o to 12° and reduced frequences of 0.1, 0.2 and 
0.3. For each case, the chordwise pressure distribution at the 
two spanwise stations, y/s = 0.48 and 0.80, between the computer 
(VSAERO-TS) and DFVLR tests are compared. 

For the present calculations, 6 spanwise strips of panels 
are used; 3 up to the mid-span station (y/s = .5) where the sweep 
starts and 3 on the outboard section. For most of the cases run, 
the number of panels chordwise is 20 (10 upper and 10 lower). 

The peak pressure coefficients for the real part varies from 
-15 to -40 as the angle of attack varies from 0° to 12°, while 
the peak pressures in the imaginary part are in the range + 3. 
The viscous effects are not taken into account in the present 
calculations. (Preliminary investigations of the influence of 
unsteady boundary layer displacement effect on the inboard re- 
gions of the rectangular tip indicated only a small viscous 
correction in the range of conditions covered here.) 

The correlation between calculated and experimental pres- 
sures at the outer station (y/s = .80) is very good, but this 
deteriorates at the inboard station (y/s = .48). The inner 
station is adjacent to the planform kink where the outboard sweep 
starts. An increase in panel density in the spanwise direction 
would probably improve the correlation here; however, this dis- 
crepancy does not appear for the unswept tip shapes. 
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Figure 4 (a) . • Comparison of Chordwise Pressure Distribution at 
y/s = 0.48 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 12°, a. = 0.704°, to = 0.1). 



Figure 4 (b) . Comparison of Chordwise Pressure Distribution at 
y/s = 0.80 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a„ = 12°, = 0.704°/ a = 0.1). 








ution 
and D: 
= 0.2 






Figure 6 (b) . Comparison of Chordwise Pressure Distribution at 
y/s = 0.80 between Computed (VSAERO-TS) and DFVLR 
Test. Swept Tip (a = 8°, a. = 0.707, w = 0.2). 




Figure 7(a). Comparison of Chordwise Pressure Distribution at 
' y/s = 0.48 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 0°, a. = 0.714°, w = 0.3). 








Figure 8(a). .Comparison of Chordw.ise Pressure Distribution at 
y/s = 0.48 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 4°, a. = 0.714°. w = 0.3). 



Figure 8 (b) . Comparison of Chordwise Pressure Distribution at 
y/s = 0.80 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 4°, a- = 0.714°, w = 0.3). 



Figure 9 (a) . • Comparison of Chordwise Pressure Distribution at 
y/s = 0.48 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 8°, a. = 0.717°, to = 0.3). 



Figure 9 (b) . Comparison of Chordwise Pressure Distribution at 
. y/s = 0.80 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a^ = 8°, a. = 0.717°, w = 0.3). 



Figure 10 (a) . Comparison of Chordwise Pressure Distribution at 
• y/s =0.48 between Computed (VSAERO-TS) and DFVLR 
Test, Swept Tip (a = 12°, a. = 0.714° , w = 0.3). 



Figure 10 (b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.80 between Computed (VSAERO-TS) and DFVLR 
Test. Swept Tip (a = 12°, a. = 0.704°, to = 0.3) . 



The correlation between the computed and the experimental 
values, for the most part, is very good, demonstrating that 
VSAERO-TS is an effective computational tool for the range of 
angles of attack and reduced frequencies covered in this study. 

The correlations are slightly better for higher frequency 
cases. The explanation for this is that in all these cases the 
same number of time steps (120 for 1 1/2 cycles) are used* For a 
proper comparison to the same level of computational accuracy, 
the number of time steps should be inversely proportional to the 
reduced frequency. This would maintain an essentially common 
step length size for the transport distance for wake elements. 


3-0 RECTANGULAR TIP 


A systematic parametric and sensitivity study has been con- 
ducted based on the rectangular tip case and the results are 
presented in a separate theory document (5). In the present 
document only one case is presented for this tip as the quality 
of the comparisons with experimental data is basically the same 
as for the swept tip, apart from an improvement in the semispan 
r<egion. 

Figure 11 shows a comparison between calculated (VSAERO-TS) 
amd DFVLR results for the case, a 0 = 120, ai = 2.0660, w = 0.30. 
®he chordwise pressure distributions for real and imaginary parts 
acre shown in Figure 11(a) through 11(d) at two spanwise stations, 
©.25 semispan and 0.70 semispan. The inboard station results are 
representative for all four tip shapes. The correlations with 
fcihe experimental data are generally very good in both the real 
ajnd imaginary parts. There is, however, a tendency for the 
calculated imaginary pressure distributions to remain "opened 
cut" towards the trailing edge (this is exaggerated in the 
present case due to the Cp scale being larger than the rest of 
fcfre imaginary pressure plots). This tendency, which became more 
pronounced in the tip regions of the taper and ogee planforms, is 
probably due to an incompatibility between the size of the shed- 
ding panel and the transport distance during one time step, V6t, 
points on the wake. This is discussed further in the theory 
document (5). 
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4.0 TAPERED TTP 


VSAERO-H was used to compute several cases for the tapered 
tip wing. The following cases of comparison between the computed 
and DFVLR experimental pressure distributions are presented. 


Figure No. 

a o 

(degrees) 

a i 

(degrees) 

0) 

12 

4 

0.713 

0.1 

13 

12 

0.711 

0.1 

14 

8 

0.703 

0.2 

15 

12 

0.703 

0.2 

16 

4 

0.705 

0.3 

17 

12 

0.701 

0.3 


Table 2 . Comparison of Computed and Experimental Pressure 
Distributions for the Tapered Tip Wing. 


The cases included in this section cover the angle of attack 
(steady) range of 0 to 120 and reduced frequencies of 0.1/ 0.2 
and 0.3. VSAERO-H was used to compute the theoretical results as 
VSAERO— TS required more than 120 time steps for each cycle to 
obtain a convergent solution in the tip region. This is due to 
the sensitivity of the shedding model to the relative time-step 
size size/shedding panel size mentioned in the previous section. 
This problem is discussed further in the theory document. 

For each case the computed chordwise pressure distributions 
are compared with the DFVLR data at three spanwise stations/ y/s 
= 0.25/ 0.80 and 0.95. The comparison between the theory and 
experiment is good in every case except at the tip section. 
Since the detailed modeling of the tip vortex is not included in 
VSAERO-H/ the discrepancy is an expected one. Although the 
VSAERO-TS and the DFVLR test results are not compared here for 
the tapered tip, the limited number of cases investigated using 
an impractically large number of time steps 0200) resulted in a 
fairly good comparison between theory and experiment. 
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Figure 12(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a Q = 4°, = 0.713°, w = 0.1). 
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Figure 12(b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.80 between Computed (VSAERO-H) and DFVLR 
• Test, Tapered Tip (a = 4°, a. = 0.713°, <m - 0.1 




Figure 13 (a) . 


Comparison of Chordwise Pressure Distribution at 
y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a_ = 12°, a. = 0.711°, w = 0. 


Figure 13 (b) . 


Comparison of Chordwise Pressure Distribution at 
y/s = 0.80 between Computed (VSAERO-H) and DFVLR 


Test, Tapered Tip (a = 12' 


a . = 0 .711° , w = 0 . 


13 (c) . Comparison of Chordwise Pressure Distribution at 
. y/s = 0.95 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 12°, a, = 0.711°. to = 0. 


Figure 14(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a =8°, a. = 0.703°, w = 0.2 



Figure 14 (b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.80 between Computer (VSAERO-H) and DFVLR 
Test, Tapered Tip (ot =8°, a. = 0.703°, w = 0.2 



Figure 14 (c) . 


Comparison of Chordwise Pressure Distribution at 
y/s =0.95 between Computed (VSAERO-H) and DFVLR 


Test, Tapered Tip 


<a Q = 


U> 



Figure 15(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 12°, a. = 0.703°, w = 0. 




Figure 15(c). Comparison of Chordwise Pressure Distribution at 

y/s = 0.95 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 12°, a. = 0.703°, oj = 0. 



Figure 16(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 4°, a. = 0.705°, w = 0.3 




Figure 16(c). Comparison of Chordwise Pressure Distribution at 

y/s =0.95 between Computed (VSAERO-H) and DFVLR 
Test. Tapered Tip (a_ = 4°, a. = 0.705°, w = 0.3 



Figure 17(a) t . 


Comparison of Chordwise Pressure Distribution at 
y/s = 0.25 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 12°, a.. = 0.701°, to = 0 . 



Figure 17 (b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.80 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (or = 12°, a. = 0.701°, to = 0. 




Figure 17(c). Comparison of Chordwise Pressure Distribution at 

y/s = 0.95 between Computed (VSAERO-H) and DFVLR 
Test, Tapered Tip (a = 12°, a. = 0.701°, w = 0. 



5.'0 OGEE TIP 


VSAERO-H was used to compute several cases for an ogee tip 
wing. The following cases of comparison between the computed and 
the DFVLR experimental pressure distribution are presented: 


gure No. 

a o 

(degrees) 

a i 

(degrees) 

to 

18 

0 

0.7 07 

0.2 

19 

12 

0.717 

0.2 

20 

4 

0.710 

0.3 

21 

12 

0.710 

0.3 


Table 3. Comparison between Computed and DFVLR Experimental 
Pressure Distribution for the Ogee Tip. 


The cases included in this section cover the range of 0 to 
12o angle of attack steady) and the reduced frequencies of 0.2 
and 0.3. For each case the chordwise pressure distribution at 
three spanwise stations, y/s = 0..39, 0.85 and 0.99, for computed 
(VSAERO-H) and DFVLR test results is compared. In all cases 
investigated, the comparison is fair. 
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Figure 18 (a) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.39 between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 0°, a. = 0.707°, a = 0.2). 



Figure 18 (b) . CompaarJison of Chordwise Pressure Distribution at 

y/ s =-. ©.85 between Computed (VSAERO-H) and DFVLR 
Test,, (Ogee Tip (a = 0°, a. = 0.707°, w = 0.2). 



Figure 18 (c) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.99 between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 0°, a. = 0.707°, a> = 0.3). 



Figure 19(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.39 between Computed (VSAERO-H) and DFVLR 
'Test, Ogee Tip (a = 12°, a. = 0.717°, to = 0.2). 



Figure 19 (b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.85 between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 12°, a. = 0.717°, w = 0.2). 




Figure 20(a). Qnnnparison of Chordwise Pressure Distribution at 
. y/ss = 0.39 between Computed (VSAERO-H) and DFVLR 
• Tesst, Ogee Tip (a = 4°, a. = 0.710°, w = 0.3). 




Figure 20 (b) . Comparison of Chordwise Pressure Distribution at 

y/s = 0.85. between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 4°, a. = 0.710°, w = 0.3). 





Figure 21(a). Comparison of Chordwise Pressure Distribution at 

y/s = 0.39 between Computed' (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 12°, a. = 0.710°, w = 0.3). 



Figure 20 (b) . Comparison of Chordwise Pressure Distribution at 
• y/s = 0.85 between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 12°, a. = 0.710°, w = 0.3). 




Figure 20(c) . 


Comparison of Chordwise Pressure Distribution at 
y/s = 0.99 between Computed (VSAERO-H) and DFVLR 
Test, Ogee Tip (a = 12°, a.. = 0.710°, to = 0.3). 


6.0 CONCLUSIONS 


A comparison between the computed (VSAERO-TS and VSAERO-H) 
and DFVLR test results for chordwise pressure distributions for 
rectangular, swept, taper and ogee blade tips is presented in 
this report. A complete discussion on the theory, limitations 
and the convergence characteristics of the VSAERO-TS and VSAERO-H 
codes are presented in a separate theory document. 

A wide range of angles of attack (mean) from 0 to 12 degrees 
and reduced frequencies of 0.1, 0.2 and 0.3 are covered in this 
report. Also, the comparison includes several spanwise stations. 
For the most part the comparison between the theory and experi- 
ment is very good for the range of conditions covered; however, 
the time-stepping calculations (VSAERO-TS) showed a sensitivity 
to the relative wake-shedding panel size and time step size. 
Results from this program were not, therefore, presented for the 
tapered and ogee tip cases. Although some good comparisons were 
obtained for these cases, the number of time steps required was 
impr actically large. Ways of alleviating this sensitivity are 
being examined. The harmonic wake calculations (VSAERO-H) were 
good over the full range of conditions considered. 
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